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DEVELOPMENT  OF  FINITE  DYNAMIC  ELEMENTS 


Two  new  dynamic  elements  have  been  successfully  developed  during  the  grant 
period.  Thus,  a  rectangular  plane  stress/strain  element  was  firstly  developed 
to  model  two-dimensional  elements.  Next  a  general  triangular  plane 
stress/strain  element  was  also  generated  to  model  more  complex  structural 
geometry. 

The  two  elements  were  extensively  used  to  model  a  wide  variety  of  structures 
to  afford  a  comparison  of  relative  performances  with  the  usual  finite  element 
method.  Such  results  clearly  indicate  that  dynamic  elements  are  much  more 
efficient  than  the  commonly  used  finite  elements.  The  two  major  benefits  of 
the  dynamic  element  method  are  summarized  below: 

1.  For  the  same  solution  accuracy  a  dynamic  element  modeling  effects  a 
considerable  reduction  in  the  number  of  structural  elements  compared 
to  the  finite  element  model,  thereby  reducing  the  data  preparation 
effort,  considerably. 

2.  Most  importantly,  the  dynamic  element  method  is  significantly  more 
efficient  than  the  finite  element  method,  effecting  a  saving  in 
solution  effort  by  a  factor  between  2  and  5. 

Results  of  these  research  efforts  have  been  published  in  articles  1,  2,4  and  5 
of  the  list  provided  later.  These  publications  include  complete  FORTRAN  IV 
program  listing  of  the  dynamic  elements  for  ready  use. 


DEVELOPMENT  OF  A  GENERALIZED  EIGENPROBLEM  SOLUTION  ROUTINE 

A  new  unified  solution  algorithm  and  an  associated  computer  program  has  been 
developed  for  the  efficient  and  economical  solution  of  a  variety  of  structural 
eigenvalue  problems,  including  the  quadratic  matrix  equations  associated  with 
a  dynamic  element  formulation.  The  program  is  capable  of  solving  a  variety  of 
problems  as  follows: 

1.  Free  vibration  analysis  of  undamped  structures. 

2.  Buckling  analysis  of  structures. 

3.  Free  vibration  analysis  of  structures  with  viscous  and/or  structural 
damping. 

4.  Solution  of  quadratic  matrix  equations  for  free  vibration  analysis 
employing  a  dynamic  element  modeling. 

5.  Free  vibration  analysis  of  undamped  spinning  structures. 

t 

6.  Free  vibration  analysis  of  spinning  structures  with  viscous  and/or 
structural  damping. 

7.  Solution  of  simultaneous  equations,  real  and  hermitian  matrices. 
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The  program  EIGSOL  based  on  a  combined  Sturm  sequence  and  inverse  iteration 
method  automatically  switches  on  to  either  the  real  or  complex  mode  of 
operation,  depending  on  the  problem  type  that  a  user  intends  to  solve.  Other 
pertinent  features  of  the  program  are  summarized  as  below: 

i.  The  program  fully  exploits  the  matrix  sparsity  such  as 

usually  associated  with  a  finite  and  dynamic  element  modeling, 
requiring  only  the  upper  symmetric  halves  of  the  relevant  matrices. 

11.  An  out-of-core  solution  strategy,  involving  only  a  moderate  amount  of 
core  storage,  enables  solution  of  complex  practical  problems  of 
almost  any  magnitude. 

iii.  A  special  symmetric  matrix  decomposition  scheme  has  been  utilized  by 
the  program  that  effects  rather  fast  convergence  to  the  roots  and 
vectors.  Also,  the  algorithm  makes  special  provisions  to  ensure 
nunerical  solution  accuracy. 

1v.  The  program  enables  computation  of  a  few  desired  roots  and  vectors 
only,  without  having  to  compute  any  others. 

v.  A  novel  bisection  scheme  has  also  been  originated  that  effects  much 
faster  convergence  to  multiple  roots,  as  compared  to  the  usual 
techniques. 

Details  of  this  new  algorithm  and  the  program  EIGSOL  are  published  in  articles 
3  and  5  of  the  list,  presented  next. 


PUBLICATIONS 

The  following  publications  provide  the  results  of  research  effort  conducted 
over  three  years  In  connection  with  the  AFOSR  grant. 

1.  "Development  of  a  Finite  Dynamic  Element  for  Free  Vibration  Analysis 
of  Two-Dimensional  Structures",  International  Journal  for  Numerical 
Methods  in  Engineering,  Vol.  12,  pp  1311-1327  (1978) 

2.  "Finite  Dynamic  Element  Formulation  for  a  Plane  Triangular  Element", 

International  Journal  for  Nunerical  Methods  in  Enqineerinq.  Vol.  14, 

ppl*3PWffTT979J -  - - 

3.  "Development  of  a  Unified  Numerical  Procedure  for  Free  Vibration 

Analysis  of  Structures".Internat1onal  Journal  for  Numerical  Methods 
in  Engineering,  ,to  be  published  In  October,  1980.  ~ 

4.  "Dynamics  of  Multibody  Systems",  invited  paper.  Proceedings  of  IUTAM 
Symposium,  held  In  Munich/Germany,  August  29-September  3,  1977. 

5.  "On 'Some  Recent  Developments  In  Nunerical  Methods  In  Structural 
Dynamics”,  Invited  paper.  Proceedings  of  "U.S. -Japan  Seminar  on 
Interdisciplinary  Finite  Element  Analysis",  Cornell  University, 

August  1978. 


-2- 


CONCLUDING  REMARKS 


Ouring  the  duration  of  the  3  year  grant  period,  most  significant  advances  have 
been  made  in  the  area  of  mmerical  analysis  pertaining  to  structural 
dynamics.  New  dynamic  element  formulations  have  been  achieved  along  with  the 
development  of  a  new  elaenproblem  solution  algorithm  and  an  associated 
computer  program.  Results  of  these  research  effort  have  been  published  in 
important  open  literature.  The  associated  computer  programs,  available  for 
ready  use,  is  expected  to  be  widely  used  by  the  mathematical  and  engineering 
community. 
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